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Electrostatic interaction is the driving force for the encapsulation by virus coat proteins of 
nanoparticles such as quantum dots, gold particles and magnetic beads for, e.g., imaging and ther- 
apeutic purposes. In recent experimental work, Daniel et al. [^(75 Nano 4 (2010), 3853-3860] 
found the encapsulation efficiency to sensitively depend on the interplay between the surface charge 
density of negatively charged gold nanoparticles and the number of positive charges on the RNA 
binding domains of the proteins. Surprisingly, these experiments reveal that despite the highly 
cooperative nature of the co-assembly at low pH, the efficiency of encapsulation is a gradual func- 
tion of their surface charge density. We present a simple all-or-nothing mass action law combined 
with an electrostatic interaction model to explain the experiments. We find quantitative agreement 
with experimental observations, supporting the existence of a natural statistical charge distribution 
between nanoparticles. 

PACS numbers: 87.15. kt,87.16.dj,34.20.-b 
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I. INTRODUCTION 

In vitro reconstitution of single-stranded RNA viruses 
is driven by electrostatic interactions between positively 
charged RNA binding domains on the virus coat pro- 
teins also known as Arginine Rich Motifs or ARMs, and 
the negative charges on the backbone of RNA molecules. 
It is not entirely surprising then, that also non-native 
RNAs can be encapsulated by virus coat proteins, as 
synthetic polyanions, supramolecular polyanions and sur- 
face functionalized nanoparticles [0-0. Indeed, nega- 
tively charged gold, silica oxide, magnetite, as well as 
zinc and cadmium sulfide nanoparticles have been suc- 
cessfully encapsulated for the purpose of manufacturing 
mcta materials and for imaging and therapeutic applica- 
tions 0, |l|-|2|. 

It is important to note in this context that the coat pro- 
teins of brome mosaic virus (BMV) and cowpea chlorotic 
mottle virus (CCMV) have been shown to spontaneously 
self-assemble in vitro into empty shells (capsids) under 
conditions of sufficiently low pH and high ionic strength 
in the absence of viral RNA Under condi- 

tions of neutral or slightly basic pHs, virus-like particles 
form only in the presence of viral RNAs or other anionic 
cargos. Plausibly, electrostatic repulsion by the RNA 
binding domains, the ARMs, on the coat proteins pre- 
vents the self-assembly of empty shells, which would be 
of evolutionary advantage to the virus [0-0 |o|-^ . 

It turns out that not only the solution conditions are 
important but also the assembly protocol. This is illus- 
trated in Figure 1, showing the reconstitution phase dia- 
gram of BMV from its constituent RNA and coat proteins 
as a function of pH and ionic strength obtained by Cuillel 
et al. by means of neutron scattering experiments [^. 



Starting off from the same initial conditions, indicated in 
the figure by a star, different structures form depending 
on the assembly protocol followed. If the pH is lowered 
first at high ionic strength then empty capsids form even 
in the presence of RNA, and filled ones are not found 
upon subsequent decrease of the ionic strength. How- 
ever, if the ionic strength is lowered first and the pH is 
lowered after that, virus- like shells with RNA packaged 
inside them do form. It shows that once shells form, 
they do not easily disassemble and reassemble into a new 
structure even if thermodynamically this would be ad- 
vantageous. Hysteresis effects observed in assembly and 
disassembly experiments confirm this [|3^-|35[|. 

To obtain a better understanding of virus assembly 
in general and of the encapsulation of nanoparticles by 
virus coat proteins in particular relevant to many prac- 
tical applications already alluded to, it is of interest to 
investigate the physical principles that determine the en- 
capsulation efficiency and how this ties in with the assem- 
bly pathways. For this purpose, Daniel et al. have re- 
cently systematically investigated the impact of the sur- 
face charge density of functionalized gold nanoparticles 
on the formation of virus-like particles (VLPs) following 
the assembly protocol indicated in Fig. by the thick 
line. The encapsulation experiments have revealed the 
following [|36|-|39|: 

1) The size of the encapsulated nanoparticles directs 
the size of the virus-like particles albeit that their 
sizes do seem to correspond to the usual Caspar- 
Klug T numbers for viruses [^, ^ ; 

2) Nanoparticles only partially covered by viral coat 
proteins arc rare at low pH. By far most particles 
are either completely covered or remain naked, in- 
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FIG. 1: Phase diagram for the in vitro reconstitution of 
Brome Mosaic Virus BMV from viral RNA and coat protein, 
adapted from [^. The asterisk indicates the starting ionic 
strength and pH, and the arrows show the assembly proto- 
col employed in Ref. for the successful encapsulation of 
gold nanoparticles by BMV coat proteins. The symbols 0, 
♦ represent free coat proteins, empty capsids, and virus-like 
particles observed in solution respectively. All the measure- 
ments associated with the efficiency of the encapsulation of 
gold nanoparticles by BMV coat proteins discussed in the 
main text were performed at the point d). 

dicating a high level of co-operativity of protein 
binding at low pHs; 

3) The encapsulation efficiency, defined as the frac- 
tion of particles encapsulated to the total number 
of particles, increases monotonically with the sur- 
face charge density of the nanoparticles and with 
the ratio of protein to particle concentration; 

4) A minimal negative surface charge density is re- 
quired for encapsulation to be observed; 

5) All other things remaining equal, the largest en- 
capsulation efficiency appears to be obtained for 
nanoparticle sizes that their assemblies yield native 
virus capsid size; 

6) Under conditions of low pH, co-assembly of 
nanoparticles and coat proteins is highly co- 
operative; yet, the efficiency of encapsulation is a 
gradual function of the surface charge density of 
the nanoparticles. 

Theoretical studies by Hagan and by Siber et al. 
1^, H show that these findings can be rationalized within 
the context of the law of mass action, describing the 
co-assembly of proteins and nanoparticles, and Poisson- 
Boltzmann theory that describes the role of electrostatics 
in stabilizing the complexes. These and other works, in- 
cluding that of Prinsen et al. on the closely related 
problem of the stability of multishell capsids, show that 



surface charge density, nanoparticle size as well as the 
number of positive charges on and length of the ARMs 
dictate the encapsulation efficiency. 

While qualitatively in agreement with experimental 
data, the cited theories cannot explain the combination of 
high co-operativity of the assembly and the rather grad- 
ual increase in encapsulation efficiency with increasing 
surface charge density [Q. In Ref. [|l2| the discrep- 
ancy between theory and experiment was explained in 
terms of metastablc states and kinetic traps, i.e., issues 
relating to the assembly kinetics. Here, we show that 
this discrepancy need not be of kinetic origin, but may 
be due to the existence of a broad distribution of sur- 
face charge densities of the nanoparticles associated with 
the way these surfaces arc functionalized. Indeed, of- 
ten nanoparticles are functionalized by attaching weakly 
acidic ligands to their surfaces, either for colloidal stabil- 
ity purposes, binding to assays or further biofunctional- 
ization/compatibilization p3[ . 

To control the surface charge density, Dragnea and col- 
laborators grafted short polymer molecules, polyethylene 
glycol or PEG, with charged and uncharged end groups 
onto the surface of gold nanoparticles ^ . The charge 
density was controlled by varying the fraction of PEGs 
with acidic carboxyl end groups relative to those with 
neutral hydroxyl end groups. In fact, at fixed overall 
grafting density the surface charge density is determined 
not only by the fraction of end-grafted polymers with car- 
boxyl end groups but also by the fraction of the weakly 
acidic carboxyl groups that are actually dissociated. For 
pHs where the assembly experiments were done, near the 
apparent pK^ value of the weakly acidic surface carboxyl 
groups |45| , we expect on average about half the car- 
boxyl groups to be dissociated. Note that the apparent 
pKa depends, in particular, on the size of the nanoparti- 



cles and the concentration of mobile ions 46 



It is important to realize that the statistical variation 
of the fraction of ionizable groups around the average 
value is the largest for pHs around their pK^ value. In 
macroscopic systems this is of little consequence but since 
we arc dealing with nanoparticles, the number of car- 
boxyl groups on them is not large and hence the varia- 
tion around the mean must be important. Interestingly, 
accounting for this natural variability of surface charge 
density of the functionalized particles we find a very good 
match between theory and experiment, if we consider a 
Gaussian distribution around the average surface charge 
density of nanosphcrcs and set the maximum charge den- 
sity of nanoparticles to what is reported experimentally. 

The mass action model that we invoke to describe the 
data is similar to that of earlier work, assuming (i) that 
the (effective) line tension associated with incomplete 
capsids is pronounced enough to prevent formation of sig- 
nificant fractions of intermediate states, and (ii) that the 
proteins do not form empty shells under the conditions 
that the experiments were done Hence, we consider 
the nanoparticles can assume only two states, i.e., either 
completely covered or completely devoid of proteins, and 
calculate the ratio of encapsulated nanoparticles to the 
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total number of nanoparticlcs as a function of their mean 
charge density. 

For simphcity, we employ Debye-Hiickel theory to de- 
scribe the electrostatic interactions between the various 
species, and for the width of the charge distribution of 
the nanoparticlcs we use a simple Langmuir adsorption 
model. Obviously, both non-linear Coulomb and charge 
regulation affect the numerical values in our calculations 
P6| . However, they turn out not change the principle 
that we focus attention on, which is the natural fluctu- 
ation in charge density of nanoparticlcs We find 
that while the stoichiometry ratio of nanoparticle and 
protein concentrations is an important factor that may 
make the all-or-nothing assembly seem less co-operative 
than it really is, within our equilibrium model the ob- 
served gradual variation of the encapsulation efficiency 
with increasing mean charge density of the nanoparticlcs 
can only be explained by the natural polydispersity of 
the charge density. Indeed, the size polydispersity of the 
particles (< 10%) is too small to account for it |Q. 

It is important to emphasize that, as far as we know, 
the effect of charge variation has been previously ig- 
nored in the literature. The importance of this effect 
in the assembly of virus-like particles is due to the size of 
the nanospheres, which are much smaller than colloidal 
particles usually investigated. The existence of such a 
distribution is not all that significant for the usual col- 
loidal particles that are tens to hundreds of nanometers in 
size. It does become important if the particles are small 
enough, i.e., are on the nanometer scale, and if charge 
regulation takes place. Since surface functionalization of 
nanoparticlcs with weakly acidic ligands is common, our 
findings should be of relevance not only to the encapsu- 
lation of gold particles, but also apply to quantum dots, 
magnetic beads and so on. In general, our findings should 
carry over to the electrostatic binding of nanoparticlcs 
onto macroscopic surfaces . 

We note that for the purpose of this paper, we have 
considered that the solution of capsid proteins and cargo 
is in equilibrium and follow the law of mass action. Re- 
cently, self-assembly studies of Zlotnick have confirmed 
that the formation of a number of viruses follow a re- 
versible path, suggesting that our assumption of equilib- 
rium conditions is good. While the last step of capsid 
formation could be irreversible based on the studies of 
Zlotnick in Ref. [|o), at the time scales relevant to the 
experiments, Zlotnick has shown that the concentration 
of proteins subunits and capsids approximately follow the 
law of mass action despite the possibility of irreversible 
steps towards the completion of full capsids. 

The remainder of this paper is organized as follows. In 
Section ||, we present our "two-state" assembly model 
and derive the relevant formulas to calculate the effi- 
ciency of formation of virus-like particles in mixtures of 
nanoparticlcs and virus coat proteins as a function of var- 
ious free binding energies and stoichiometry ratios. The 
effect of a fixed and a Gaussian surface charge distribu- 
tions of the nanospheres is demonstrated in Sections [II 
and [V, respectively. Numerical results are illustrated 



m Section We summarize our work and present our 
main conclusions in Section VL where also discuss in 



more detail the justification for employing the approx- 
imate Debye-Hiickel theory. Details of our calculations 
are relegated to the appendix. 



II. ENCAPSULATION OF NANOSPHERES 

As explained in the introduction, the experiments de- 
scribed in Ref. |Q were performed under the conditions 
where virus-like particles do not form if their charge den- 
sity is below a certain critical charge density-see item 4 
in the previous section. This indicates that in these ex- 
periments, the concentration of protein subunits, c, must 
be below the critical subunit concentration, c,, the mini- 
mum concentration necessary for formation of empty cap- 
sids without polyanionic cargos. Otherwise, one would 
expect that the filled or empty capsids appear in solu- 
tions regardless of charge density of nanoparticlcs. Fur- 
ther, the experiments reveal that the assembly process at 
pH 4.5 is highly cooperative: the nanoparticlcs are either 
naked or fully covered by protein subunits |3^ . 

Based on these experimental observations, we consider 
a dilute solution of capsid proteins with concentration c, 
and M negatively charged nanospheres. To mimic the 
experimental conditions, we assume that c < c* and thus 
there are no empty capsid in solution. Further, we con- 
sider that nanoparticlcs can assume only two states: ei- 
ther they are fully covered with q protein subunits or 
are completely naked. Let /S.g < denote the binding 
free energy of a protein subunit adsorbed on a covered 
nanosphere. It accounts for the attractive hydrophobic 
interaction that drives the assembly in the absence of 
payload, the repulsive electrostatic interactions between 
proteins as well as the attractive electrostatic interaction 
between the proteins and the charged nanospheres. 

The partition function of a semi-grand ensemble of 
coat proteins on the collection of M empty and covered 
nanoparticlcs can then be written as 



Z = 



M 

E 

Af=0 



Ml 



- (1 + A) 



Nl{M - N)'. 



M 



(1) 



with A = cxp [/3{fi — Ag)q] the fugacity of the coat pro- 
teins on the nanoparticlcs, 13 = l/ksT the reciprocal 
thermal energy, fc^ Boltzmann's constant and T the ab- 
solute temperature. Furthermore, q is the number of pro- 
teins on a virus-like particle (VLP) and fi the chemical 
potential of the proteins. The VLPS seem to obey the 
Caspar-Klug T-numbers, and thus we have g = 60 x T 
pO| , |4l[ | . Assuming the chemical equilibrium between the 
proteins in the VLPs and those in free solution, the chem- 
ical potentials must be equal. 

The binomial distribution accounts for the number of 
indistinguishable ways of covering iV = 0, 1, M out of 
A/ nanoparticlcs. How many of the M particles arc on 
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average covered by protein depends on the solution con- 
ditions through the chemical potential ^ and binding free 
energy Ag. The efficiency rj defined as the ratio of the 
number of nanospheres encapsulated by protein subunits 
to the total number of nanospheres is now straightfor- 
wardly calculated, 



A 



(N) _ d\nZ 
~W ~ ~ ainA ~ 1 + A' 



(2) 



where {N) is the thermal average of the number of fully 
formed capsids. This is the familiar Langmuir adsorption 
isotherm. 

In view of our discussion above, a quantity of partic- 
ular interest is the strength of attractive interaction be- 
tween the negatively charged nanospheres and protein 
subunits. It depends directly on the surface charge den- 
sity of nanoparticles and drives the co-assembly. In the 
experiments of Rcf. , the nanospheres (the gold par- 
ticles) are coated with two different types of ligand: end 
carboxylated PEG (PEG-COOH) and end hydroxylated 
PEG (PEG-OH). The latter is neutral under aU condi- 
tions of pH, while the former is a weak acid and hence can 
acquire a negative charge. For a given pH, the carboxyl 
group COOH can exist in two states of protonation: pro- 
tonated and hence electrically neutral or deprotonated 
and negatively charged, depending on the acid constant 
or the negative logarithm of it, the pKa- Hence, the sur- 
face charge density of nanospheres strongly depends on 
the relative numbers of the two ligands per unit area. 

In the following section, we first assume that all 
nanospheres have the same surface charge density that 
depends to some extent on the pH of the solution. In the 
section following that, we consider the case in which the 
surface charge density varies among particles and obey a 
Gaussian distribution around an average value dictated 
again by the pH . The width of this distribution is calcu- 
lated from a simple association-dissociation model. 



III. FIXED CHARGE DENSITY 

To calculate the efficiency of encapsulated nanoparti- 
cles, it is advantageous to rewrite Eq. ^ as a function of 
two relevant binding constants: 



and 



= exp(-/3A5) 



-f = exp (-13 Ago) 



(3) 



(4) 



where Ago is the binding free energy associated with a 
subunit in an empty capsid. The latter describes the ef- 
fect of attractive hydrophobic and repulsive electrostatic 
interactions between the coat proteins in an empty shells. 
We presume these interactions to be the same in empty 
and filled shells. Note that the electrostatic repulsion is 
dominated by the positive charges on the ARMs pl| . 




FIG. 2: Illustration of a virus-like-particle. The surface 
of inner shells is negatively charged. The N-terminal tails 
(the RNA binding domains) are uniformly distributed in the 
darker shell and are positively charged. 

The binding free energy difference between filled and 
empty shells can now be written as 



qAg - qAgo = A/e, 



which gives 



K 



cxp(-/3A/s), 



(5) 



(6) 



with /e the electrostatic attractive interaction between 
a nanosphere and a protein subunit. 

The attractive electrostatic interaction between 
nanoparticles and coat proteins is the result of negative 
charges on nanospheres and the positive charges on ARM 
regions of proteins. As explained in the introduction, the 
ARM region is a polypeptide chain that extends into the 
capsid interior and is rich in arginine and lysine, amino 
acids with strongly basic residues. 

To calculate the electrostatic interaction between the 
protein shell and nanoparticles, we for simplicity assume 
that the charges on the ARMs arc uniformly distributed 
in the space between the inner surface of the protein shell 
and nanoparticles, see Fig. |2| We then solve the relevant 
linearized Poisson-Boltzmann equation and find the elec- 
trostatic energy to obey 



167r^ABO'i/02^i ((t-Ri 



+ l)_gK(fll-i?^2)(^^2 + l)) 



«;2 {kRi + 1) 



(7) 



with CTi < the nanosphere surface charge density, i?i 
the nanosphere radius plus the length of the ligands, i?2 
the inner radius of capsid and p2 > its effective volume 
charge density due to the presence of the ARMs. The 
Bjerrum length, Xg = e'^ /{AirekBT), where ethe dielec- 
tric permittivity of the medium, is 0.7 nm for water at 
room temperature, k = ^JEttXdPs is the reciprocal of the 
Debye screening length with ps the number density of 1- 
1 electrolyte. At room temperature it is about 0.3/^^, 
with Cs the molar concentration of the electrolyte. Equa- 
tion 1^ is derived in Appendix 

Further headway can be made by noting that the re- 
ciprocal of the binding free energy /Cq is related to the 
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critical capsid concentration, c,i; = with v the 

"interaction" volume, that we take to be the volume of a 
solvent molecule making c,w a mole fraction. The critical 
capsid density c* is a very important physical quantity 
that can be measured experimentally [^T| , |5^ . Equation 
can then be written as 



exp(^Ag) = 



K 



c^vKq 



(8) 



where the quantity K/Kq can be obtained from Eqs. 
and0). 

The expression for the efficiency, Eq. ^ can be simpli- 
fied further if we consider that the system is in chemical 
equilibrium and that the solution is dilute. The chemical 
potential /x of the free subunits is 



/3/i = In c/w, 



(9) 



with Cf the concentration of free subunits, which by mass 
conservation must be equal to 



qriCp = c(l - rir), 



(10) 



where c is the overall concentration of protein subunits 
and Cp that of the nanospheres in solution; the stoichiom- 
etry ratio r is defined as qcp/c so that for a complete VLP 
r = 1. Note that here we ignore a reference chemical po- 
tential. This merely renormalises Kq and is implicit in 
the value of c, . 

Substituting Eqs. (|), (|) and (0) into Eq. (|), we find 
the following expression for the efficiency rj, 



(11) 



This is our main result corresponding to the encapsu- 
lation efficiency of nanoparticles with a given surface 
charge density. 

In the limiting case K/Kq — > 1 and c < c*, we find 
that the efficiency goes to zero, meaning that weak elec- 
trostatic interaction between the nanoparticles and pro- 
tein subunits results into low efficiency of the formation 
of VLPs, as observed in the experiments [Q. For large 
enough values of K/Kq > c*/c > 1, corresponding to 
strong attraction between nanospheres and coat protein, 
the efRciency goes toward unity as expected. In Fig. ^, 
we plot the efficiency as a function of protein concen- 
tration scaled to the critical capsid concentration, c/c* 
for varying surface charge densities ai and stoichiometry 
ratios r. Regardless of the values of r the electrostatic 
interactions determine where the efficiency begins to rise 
from 0. As illustrated in the figure, the efficiency of VLPs 
gets smaller, as stoichiometry ratio becomes goes larger. 



IV. VARIABLE CHARGE DENSITY 

To include in the theory the charge density distribu- 
tion of nanoparticles, we explicitly consider the acid-base 



equilibrium of the carboxyl groups attached via the PEGs 
to the surface of the nanoparticles. The reaction scheme 
for a carboxylated PEG can be written as 

PEG-COOH ^ PEG-COO-+H+ 

The (grand canonical) partition function Zh associated 
with this equilibrium reads 



Zh = 1 + X 



H 



(12) 



with Xh = exp^PfiH — P^h) the fugacity of the hydrogen 
ions in solution, their chemical potential and — > 
the free energy of ionization. Note that because of chem- 
ical equilibrium, the chemical potential of the hydrogens 
on the PEG-COOHs and those in solution must be equal 
- the ones in solution act as a reservoir. 

In the partition function Eq. (|l^, the unity in the 
sum arises from the deprotonated state, PEG-COO~, 
and the term \h corresponds to the protonated state, 
PEG-COOH. The thermal average of the occupancy of 
H+ binding sites PEG-COO-\ forming PEG-COOHs, 
can be calculated from 



(n) 



dlnZu 
9 In A 



A 



H 



H 



1-f A 



H 



(13) 



which is equivalent to a Langmuir adsorption isotherm. 
The variance of the occupancy can be obtained from 



A 



H 



d{n) 

dlnXH (1 + Aa)2 



(14) 



The fugacity Xh can be rewritten in terms of pH and 
pKa of the weak acid as follows [^Sj, 



A 



H 



10 



-pH+pKa 



(15) 
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FIG. 3: The encapsulation efficiency r\ versus the total protein 
concentration c, scaled to the critical capsid concentration 
c* for surface charge densities a\ = 0.05 (dashed lines) and 
0.01 (solid lines) e/nm?, with three different values of the 
stoichiometric ratio r = 0.1, 1, and 10. For these plots, we 
have set the radius of the nanoparticles to Tii = 7.9 nm, that 
of the cavity of the capsid to R2 ~ 11.4 nm, q = 90 (90 
dimers assembling around the core) and the inverse Debye 
length K = 0.89 nm"\ We assume all PEG-COOH groups 
are deprotonated and negatively charged. 
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because, apart from a standard chemical potential that 
we incorporate in the ionization free energy , l3fJ,H = 
Iu^chv) with chv the mole fraction of hydrogen ions in 
the solution. Note that pH = — logjg(55.6 x chv) and 
pKa = — logiQ(55.6 X exp(— /3e//)), where we used the 
molarity of water, 55.6 M. Substituting Eq. ( |l5|) into 
Eq. (O), wc obtain 



1 



1 + IQPH-pKa 



(16) 



with a similar expression for (Ar?). This IS m essence 
the familiar Henderson-Hasselbalch equation. 

We have to translate this expectation value for a sin- 
gle dissociable group to that for a collection of them on 
the surface of the nanoparticles. There, the effective pKa 
will depend on the strength of the electrostatic interac- 
tion between them, potential hydrogen bonding between 
neighboring COOH and C00~ groups, and so on. In 
fact, the average protonation of the carboxyl groups on 
naked particles may well be different from that of encap- 
sulated ones due to charge regulation Q, This is 
an important issue but does not qualitatively change our 
results. The focus of this paper is to show that statisti- 
cal fluctuations of the charged state of the particles has 
a large impact on the perceived level of cooperativity of 
the encapsulation process. 

For simplicity we, thus, assume that there exists a sin- 
gle pKa for both naked and covered particles that de- 
pends on the concentration of salt, the surface coverage 
of PEG-COOHs and some other factors. This implies 
that we ignore any coupling between the charged states 
of the carboxyl groups on the particles 0. Grant- 
ing this simplification, which implies that our pK^ is an 
effective one, the average surface charge density can be 
found by 



1 



(17) 



with amax the maximum surface charge density of the 
acidic groups that depends on the overall grafting den- 
sity of PEGs and the ratio of the number of PEGs with 
COOH to those with OH end groups. This quantity is 
an experimentally known quantity |Q. Note that the 
PEG-OHs do not contribute to the charge density other 
than diluting the PEG-COOHs resulting into a decrease 



The variance, then, is 



^Q—pH+pKa 



qi (1 + \Q-pH+pKay 



(18) 



where is the number of all the PEG-COOHs on the 
nanosphcre. In the experiments of |3^, the maximum 
value of is ^3000 for a nanosphcre with a diameter 
of 12 nm. In principle, wc would like to know the full 
distribution P{(j), rather than the first two moments of 
that distribution. However, because the large number 
of nanoparticles and the dilute solution, we invoke the 
central limit theorem assuming that the charge distribu- 
tion is Gaussian. Hence, the probability density to find a 



nanosphcre with a surface charge density between cr and 
a + 6a, with 6 a an infinitesimal increment, is 



P{a) 



1 



V27rA<T' 



: CXp 



2AO-2 



(19) 



The theory of the previous section for a fixed charge 
density needs to be modified to accommodate the charge 
distribution on the particles. To this end, we modify the 
conservation of mass, Eq. [10, 



qcp 



daP{(T)-q{a) 



(20) 



1] 



where rji^cr) is the efficiency of encapsulation of nano 
spheres for a given charge density cr as described by 
Eq. (prj). The upper limit of the integral in Eq. ^ is 
set to the maximum value of the surface charge density, 
'^max- The value of amax is chosen assuming that all 
the -COOH groups are negatively charged, and is used 
in Eqs. 17 and ^ to calculate the mean {{a)) and the 
variance ((Act)) charge density. Due to small values of 
the variance ((Act)), the portion of the Gaussian distri- 
bution corresponding to ct < is almost zero, and thus 
the lower limit of integral in Eq. ^ is set to zero. Sub- 
stituting Eq. (11) into Eq. (po|), we obtain 



c 



daP{a) 



(jlM 

^ c* Kq (ct) c 



(21) 



The binding constants in Eq. ( ^ ) depend on the sur- 
face charge density of nanospheres. For a given c/c* < 1 
and the stoichiometry ratio r, Eq. (|2|) can be solved self- 
consistently and the ratio c//c can be determined. The 
total efficiency of encapsulation, 77, is then obtained by 
the integral given in Eq. (^). 

In Fig. ^, wc plot the efficiency as a function of the sur- 
face charge density at pH — pKci=:4.6 and pH ^ pKa=0 
for r = 0.1, 1, and 10. The efficiency curves for pH ^ 
pKa=0 grow much faster and agree with the results for 
the fixed charge density model of the previous section, as 
expected. 



V. COMPARISON WITH EXPERIMENTS 

As noted in the introduction, the efficiency of the en- 
capsulation of functionalized nanoparticles by virus coat 
proteins was systematically measured as a function of 
their surface charge density by Dragnea and collabora- 
tors p^ . The experiments were performed at pH = 4.6 
and low ionic strength. The protein concentration and 
the assembly protocol, indicated by the thick line in Fig. 
1, were chosen such that no empty capsids were formed 
in the solution regardless of the charge density of gold 
nanoparticles. If the charge density of the nanoparticles 
was increased to a value above a critical charge density, 
capsids filled with nanospheres appeared in solution. To 
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describe these experimental observations, we employ the 
two-state model presented in the previous section. 

According to the two-state model, the protein subunits 
are either free in solution or belong to VLPs, consistent 
with the experimental findings. To compare our results 
with the experiments, we need to know the values of the 
following quantities: the inverse Debye length k, the sto- 
ichiometric ratio r, the ratio of protein concentration to 
the critical concentration c/c*, the number of protein 
subunits in each VLP, q, the radii i?i and R2, the ef- 
fective volume charge density due to the presence of the 
ARMs p2, the maximum surface charge density of the 
particles cjmax and, finally, the effective pKa of the car- 
boxyl groups and the pH of the solution. 

Unfortunately, there is not sufficient experimental data 
available to fix the values of all of these parameters. In 
particular, the values for pKa and c/c* are unknown 
for the conditions under which the experiments were 
done. Therefore, we use physically reasonable estimates 
to match our theoretical predictions with the experimen- 
tal data associated with the assembly of gold nanoparti- 
cles by BMV capsid protein subunits |^ . Note also that 
there is some freedom in the choice of R2 too, because 
the inner surface of the capsid is not actually smooth but 
quite lumpy. 

Figure || illustrates the efficiency vs. surface charge 
density obtained in the experiments (empty circles) for 
the radii i?i = 7.9 nm and R2 = 11.4 nm, the stoichiom- 
etry ratio r = 0.75, and k = 0.89 nm~^. The dashed line 
in Fig. |[ corresponds to the efficiency obtained using the 
two-state model with the charge density of nanoparticles 
fixed. While the critical concentration, c*, for the forma- 
tion of empty capsids is not known for the pH at which 
the experiments were done, the quantity c,/c associated 
with the experiments should be less than one, c*/c < 1, 
based on the fact that no capsids formed when the sur- 
face charge density of nanosphcres was zero. Note that 
while the fit is not very sensitive to the small variations 
in i?2 and c/c*, the value of pK^ has a huge impact 
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0.1 0.2 0.3 0.4 0.5 

Surface charge density <3^^^ 



FIG. 4; The encapsulation efficiency r\ versus the max- 
imum surface charge density Omax e/nm^ for two cases: 
pH=pKa=4.6 (solid lines)and pH ^ pK^ = (dashed lines) 
with three values of the stoichiometric ratio r — 0.1, 1, 10. 
The other model parameters are the same as given in Fig ^ 
Ri = 7.9 nm, R2 = 11.4 nm, c/c. = 0.14. 



on the efficiency curves. This point is illustrated in Fig. |6| 
in which several efficiency curves for different values of 
pKa are plotted. All dashed lines in the figures corre- 
spond to i?i = 7.9 nm, R2 = 11.4 nm, c/c* = 0.14 and 
all solid lines correspond to i?i = 7.2 nm, R2 ~ 11 nm, 
c/c* = 0.2. The only variable among the solid or dashed 
curves is pK^ as indicated in the figure. 

As illustrated in fig. H, the two-state model with a fixed 
charge density predicts the presence of a critical charge 
density below which no VLPs form, consistent with the 
experiments. However, it docs not capture one of the im- 
portant features of the experimental data: The efficiency 
vs. charge density curves associated with the two-state 
model rise much faster than the experimental data. To 
this end, we augment the two state model with a Gaus- 
sian charge distribution as explained in the previous sec- 
tion. 

As shown in Eq. ^ the width of Gaussian charge dis- 
tribution depends on the dissociation constant (pK^), 
which is not known for carboxyl group binding to PEG. 
Several experiments show that pK^ depends on Coulomb 
and other types of interactions between neighboring car- 
boxyl groups and also on the charge density of ARM 
regions. To fit the data, we choose pK^ = 6.3, which 
is within the reasonable range of known values for the 
pKq of carboxyl end groups Q, It is a bit below 
that of 11-mercaptoundecanoic acid end grafted to gold 
nanoparticles under similar conditions of ionic strength 
p6[ . This we attribute to the effect of charge regulation 
resulting from the interaction of the carboxylic acid with 
the ARMs in the VLPs that we do not explicitly model 
in the present paper . 

Despite the simplicity of our model, the solid lines 
in Figs. ^ show that the efficiency plots obtained with 
a Gaussian charge distribution gradually rise to unity 
rather than do so abruptly if we ignore the existence of a 
charge distribution. Indeed, there is a very good match 
between the theory and experimental data. 



Fixed-charge model — — 
Gaussian model 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Surface charge density o^^^^ 



FIG. 5: The experimental data||36|] (red dots) compared to 
the Gaussian model (green solid line) and fixed charge density 
model (blue dashed line). The fitting parameters are Ri — 7.9 
nm, R2 = 11.4 nm, c/c, = 0.14, and pKa = 6.3. The choice of 
the other model parameters is based on the experimental data: 
q = 90 (90 dimers assembling around the gold nanoparticle) , 
r = 0.75 (the stoichiometry ratio), and n = 0.89 nm~^ (the 
inverse Debye length). 
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0.5 1 1.5 2 2.5 3 3.5 
Surface charge density o^^^ 



FIG. 6: The fitting parameters for daslied lines are Ri = 
7.9 nm, R2 = 11.4 nm, c/c* — 0.14 and for solid lines are 
Ri = 7.2 nm, R2 = 11 nm, and c/c, — 0.2. The choice of the 
other model parameters is based on the experimental data: 
g = 90 (90 dimers assembling around the core), r = 0.75 (the 
stoichiometry ratio), and a = 0.89 nm~^ (the inverse Debye 
length). The value of pKa is indicated next to each curve and 
varies from 5.4 to 6.5. 



VI. SUMMARY AND CONCLUSION 

Because the structure of transient assembly interme- 
diates is not easily accessible experimentally, the under- 
lying mechanisms of virus assembly are not, in general, 
very well understood. Still, several important theoreti- 
cal and experimental works have recently shed light on 
the process of the formation in vitro of virus capsids, 
in the presence and absence of polyanionic payload 

H |0[ m, H, ||. Relevant to the work presented 
in this paper are the experimental studies of Dragnea et 
al. [|6[ ^1 associated with the efficiency of incorporation 
of nano-payloads as a function of surface charge density. 
According to their studies, there is a critical charge den- 
sity below which no encapsulated nanoparticles are ob- 
served. 

A prominent feature of the self-assembly studies in Ref. 
p6[ is the dependence of the efficiency of incorporation 
of nanoparticles on surface charge density. Dragnea et 
al. found that beyond a critical charge density, the num- 
ber of complexes of capsids and nanoparticles increases 
gradually as the charge density of nanoparticles increases. 
The gradual rise in the efficiency curves is unexpected be- 
cause the process of capsid assembly involves a high-order 
mass action that is inherently co-operative About 
180 protein subunits should assemble to form the T = 
3 structure of BMV capsids, and this is the reason why 
there is always a very sharp rise in the the efficiency of 
encapsulation of nanoparticles beyond a critical protein 
concentration. 

In order to understand the experimental results of Ref. 
p6t , Hagan developed equilibrium and kinetic theories to 
explore the role of the charged nanoparticles in the as- 
sembly of virus-like particles [|l^ . His studies nicely show 
that there is indeed a threshold surface charge density 
above which nanoparticles become efficiently encapsuled 



by viral capsid proteins. However, the theoretical studies 
of Ref. Q also show a sharp and very steep increase of 
the encapsulation efficiency with increasing charge den- 
sity. Metastable states and the presence of kinetic traps 
were suggested as part of the reason for the discrepancy 
between the theory and experiments. 

The central idea of this paper is that a natural, sta- 
tistical charge variation between different nanoparticles, 
a result of the way they are made to acquire charges 
not only in the above experiments but quite generally, 
must be the root cause of the discrepancy between previ- 
ous theory and the available experimental data . 
The effect of charge variation has, as far as we know, 
been previously ignored in the literature. We argue that 
the importance of this effect in the assembly of virus- 
like particles is due to the size of the nanospheres, which 
are much smaller than colloidal particles usually investi- 
gated. The larger the particles are, the smaller the effect 
of charge variation is, which seems to be borne out by 
the experimental data presented in Ref. Q . 

We put forward a very simple model to describe the 
assembly of nano-cargos by virus coat proteins. In par- 
ticular, we calculated the efficiency of nanospheres en- 
capsulated by capsid proteins as a function of surface 
charge density, the coat protein concentrations, the sto- 
ichiometric ratio of protein to cargo concentrations, and 
interactions between protein and cargo. Our model is 
also able to explain successfully the presence of a criti- 
cal charge density. We show that discrepancy between 
the tlrcory presented in Ref. fl^ and the experiments of 
Dragnea et al. can be explained by a distribution in the 
charge density of the nanoparticles. 

Our work shows that the existence of a broad distri- 
bution of surface charges on nanoparticles stabilized by 
weakly acidic groups, can have a strong impact on their 
encapsulation by virus coat proteins. This is particularly 
true for pHs near the pKa of such groups, where the 
charge variation is most pronounced. The experiments 
cited were indeed performed close to the expected pKa 
between 5 and 7 for end-grafted polymeric surface layers 
bearing carboxyl moieties. 

It is important to point out that our findings are not 
significantly affected by the Debye-Hiickcl approximation 
that we employed in our calculations for reasons of sim- 
plicity. Our goal was not to let complicated details of the 
theory requiring additional adjustable parameters cloud 
the issue. Indeed, including charge regulation and solving 
the non-linear Poisson-Boltzmann equation in the spirit 
of Prinsen and collaborators [Q, we find that without 
charge variation the experimental data cannot be ex- 
plained As we alluded to in the introduction, also 
more elaborate theories [|l^ that in addition model the 
conformational statistics of the PEGs cannot explain the 
gradual rise of the encapsulation efficiency with surface 
charge density. 

In view of the rapidly increasing number of self- 
assembly studies involving viral proteins and functional- 
izcd nanoparticles, we feel it is important to convey the 
message that there is an additional complication that is 
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not considered: fluctuations become more important the 
smaller the particles. Clearly, comprehensive experimen- 
tal and theoretical investigations of the parameters that 
affect capsid assembly in general, and the assembly of 
VLPs with poly- or macroanionic particles in particular, 
are highly desirable, and should have a significant im- 
pact on the development of optical imaging assays [ p2[ , 
therapies and treatment of viral infection [^8| . 

Finally we emphasize that the calculation done in this 
paper is quite general and could be applied to all nano- 
particles functionalized with weakly acidic groups. Note 
that weakly acidic groups are commonly used in the con- 
text of biocompatibalizing nanoparticles . 

The authors gratefully acknowledge helpful discussions 
with Mehran Kardar. We are grateful to Bogdan Drag- 
nea and Irena Tsvetkova for discussions and providing the 
experimental data. This work was partially supported 
by the National Science Foundation through grant No. 
DMR-06-45668. 



Appendix A: Electrostatic interaction energy 



w ith the invers e of Debye length defined as k = 
^j2NAe^I /ksTe, I the ionic strength (M), Na Avo- 
gadro's number and e the dielectric permittivity of the 
medium. 

The general solution of eletrostatic potential can be 
found to be 



_ ^Gxp_(fcr) 



B 



cxp ( — Kr) 



jj Cxp(-Kr) 



< r < i?i 
-f- P2/K^e, Ri <r < R2 
r>R2 

(A2) 

The four unknown parameters A, B, C, D can be obtained 
applying the boundary conditions at Ri and i?2- The 
resulting electrostatic free energy can then be written as 



(A3) 



To calculate the electrostatic interaction between a 
nanosphere and N-terminal tail of capsid proteins, we 
consider an electrolyte solution with dielectric constant 
e. The electrostatic potential obeys the following lin- 
earized Poisson-Boltzmann equations in three different 
regions (see Fig. 2), 



-cri(5(r - i?i)/e, 0<r<Ri 

- P2/e, Ri <r < R2 

K^(t), r > i?2 



which yields 



47re2crip2i?? {{uRi + 1) - e'^'^^^-^'^^ {nR2 + 1)) 



en^ {kRi + 1) 



(Al) and self-energy to be 



(A4) 



fself 



2TTe^alRl ne^pl ((^i?i + 1) {2k^RI + ^^^Rj (3 - 2KR2) - 3) - 3e^"(«^-^^) {kRi - 1) {kR2 + 1) - 
e{KRi + 1) 3^56 (Ki?i + 1) 



(A5) 



r 



Note that the contribution of the entropy of the ions 
is negligible within the Debye-Hiickel approximation. In 



Eq. |A4| , if CTi and p2 have opposite signs then the binding 
constant {K/KqY = exp(-/3/B) >> 1. 
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